arXivil506.00594v3 [cond-mat.soft] 5 Oct 2015 


Nucleation of symmetric domains in the coupled leaflets of a bilayer 

J. J. Williamsoni'El and P. D. Olmstedi’S 

^Department of Physies, Institute for Soft Matter Synthesis and Metrology, 

Georgetown University, 37th and O Streets, N.W., Washington, D.C. 20057, USA 

(Dated: October 6, 2015) 

We study the kinetics governing the attainment of inter-leaflet domain symmetry in a phase- 
separating amphiphilic bilayer. “Indirect” inter-leaflet coupling via hydrophobic mismatch can in¬ 
duce an instability towards a metastable pattern of locally asymmetric domains upon quenching 
from high temperature. This necessitates a nucleation step to form the conventional symmetric 
pattern of domains, which are favoured by a “direct” inter-leaflet coupling. We model the energetics 
for a symmetric domain to nucleate from the metastable state, and find that an interplay between 
hydrophobic mismatch and thickness stretching/compression causes the effective hydrophobic mis¬ 
match, and thus line tension, to depend on domain size. This leads to strong departure from classical 
nucleation theory. We speculate on implications for cell membrane rafts or clusters, whose size may 
be of similar magnitude to estimated critical radii for domain symmetry. 


I. INTRODUCTION 

A phase-separating lipid (or other amphiphilic) bilayer 
may access competing equilibrium and metastable phase 
coexistences, due to the presence of two leaflets subject 
to competing inter-leaflet couplings (Figs. [U [2^) [J, 0- 
A “direct” inter-leaflet coupling [sl-Q promotes registered 
(R) bilayer phases, in which both leaflets are locally dom¬ 
inated by the same species. An “indirect” coupling from 
hydrophobic tail length mismatch favours uniform bilayer 
thickness and thus promotes antiregistered (AR) 

phases, in which the leaflets are dominated by differ¬ 
ent species so that the bilayer is locally compositionally 
asymmetric. Depending on the choice of overall leaflet 
compositions, two (e.g., R-R or AR-AR) or three bilayer 
phases may coexist. For example, asymmetric overall 
leaflet compositions can lead to two approximately sym¬ 
metric bilayer phases and a highly as;j^metric one (R-R- 
AR), first observed and explained in j6|. Phase equilibria 
of coupled bilayer leaflets were subsecmently studied us¬ 
ing phenomenological free energies [3, S, [l2- However, a 
full description of the indirect coupling described above 
requires a model that microscopically incorporates hy¬ 
drophobic mismatch [l|. 

Using such a model, we have shown how coexistence 
of antiregistered phases can be kinetically preferred due 
to the effect of hydrophobic mismatch [l[ so that, before 
equilibrating, a quenched bilayer must escape a (typically 
metastable) locally asymmetric state. The understand¬ 
ing of this novel statistical thermodynamics will allow 
greater control over artificial membranes, and the under¬ 
lying interactions are expected to play a role in trans¬ 
membrane organisation of rafts or clusters in vivo, with 
possible relevance to signalling and anaesthetic ac¬ 
tion [T^. The predicted behaviour constitutes an exam¬ 
ple of Ostwald’s “rule of stages” M , by which a system 
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FIG. 1. Partial phase diagram showing competing equilib¬ 
rium R-R (black) and metastable AR-AR (red) phase coex¬ 
istences [l|. Spinodals enclose the regions of local stability. 
Cartoons of the dominant inter-leaflet arrangement in each 
bilayer phase are shown. The grey dotted line illustrates R- 
AR coexistence, which is briefly discussed in Section IIII Al 
Other phase coexistences not considered here are omitted 
[l[. Parameters: Ao — 2a, n — 3a~^kBT, V = 0.6/cbT, 
J = 4: a B = 0.48 a ^k^T. 


will pass through available metastable states on its way 
to equilibrium. Ostwald’s rule is familiar (via different 
origins) in colloids, metallurgy and drug design. 

Immediately after quenching a bilayer to a phase- 
separating region of parameter space, any spinodal insta¬ 
bilities to which the uniform state is subject compete to 
determine the dominant initial demixing mode. A bilayer 
with roughly equimolar composition in each leaflet [l5| 
can be subject both to an “R mode” with composition 
perturbations locally symmetric between leaflets, and a 
perpendicular AR mode with asymmetric perturbations 
(Fig. [1]). If the AR mode is fastest-growing, spinodal de- 
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composition to AR-AR coexistence occurs first, leading 
to local asymmetry throughout the bilayer. To equili¬ 
brate to R-R from a metastable AR-AR state, the bilayer 
must undergo nucleation of registered bilayer phases. 
Hence, three classes of kinetics arise: direct separation 
into equilibrium phases, equilibration via nucleation out 
of a metastable state, or trapping in a metastable state. 
For other overall compositions subject to competing in¬ 
stabilities, the competition of symmetric and asymmetric 
phases is qualitatively similar 0 , though more complex. 

In this paper, we focus on the nucleation energetics 
that determine whether equilibrium domain symmetry 
is reached from a metastable state. First, we introduce 
the model and discuss the interpretation of bilayer do¬ 
main symmetry and asymmetry via phase diagrams with 
a composition axis for each leaflet, with reference to ex¬ 
isting experiment and theory. We then identify the three 
classes of kinetics in simulation, guided by a linear in¬ 
stability analysis, and develop a theory for the nucle¬ 
ation of registered domains, which captures the interplay 
of bulk free energy with thickness mismatch occurring 
at the perimeter of a registered domain. Together with 
the linear stability analysis, the calculated nucleation en¬ 
ergetics are consistent with the simulation results. We 
find that the effective hydrophobic mismatch between a 
domain and its surroundings is domain size-dependent, 
which causes strong departure from classical nucleation 
theory. 


II. MODEL 

A detailed description of the model appears in 
We briefly recapitulate the model and its analysis in 
terms of phase diagrams and kinetics. A schematic of 
the lattice model and how it is coarse-grained is shown 
in Fig. [21 


A. Lattice model 

A local bilayer patch is described as N lattice sites 
with top (t) and bottom (b)-leaflet lipids (Fig. | 2 ] 3 ). The 
lipids’ hydrophobic lengths lead to the total bilayer 
thickness di = and inter-leaflet difference = 

Model species S and U represent saturated and 
unsaturated lipids or, e.g., the liquid-ordered (To) and 
liquid-disordered (Ld) states in a ternary mixture Q- 
The Hamiltonian is 

<i,j> <i,j> 

+ + E F ((^1 - ^0 , 

i i 

( 1 ) 

where = 1 if the top (bottom) of lattice site i con¬ 
tains an S lipid, = 0 if U. The species-dependent 
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FIG. 2. (a) Mixed bilayer containing S and U {Lo or gel 

and Ld-like) model species, illustrating the locally symmet¬ 
ric (R-R) and locally asymmetric (AR-AR) phase coexis¬ 
tences considered here, (b) Microscopic lattice model for cou¬ 
pled leaflets, which can be coarse-grained (CG) (c) to give 
the mean-field free-energy density /(0^,0^) as a function of 
locally-averaged leaflet compositions, and analysed for kinet¬ 
ics of domain formation with the inclusion of gradient costs 
for domain boundaries [i[. The lattice model can also be di¬ 
rectly simulated. 


ideal (i.e., preferred) hydrophobic tail lengths are = 
iso for an S lipid at the top (bottom) of site i, or iuo for 
I/, and each site is pairwise registered (R, SS or UU) or 
antiregistered (AR, SU or US). 

V = Vio — ^(Vbo + khi) quantifies purely intra-leaflet 
interactions, such as those between headgroups. The “di¬ 
rect” coupling B promotes pairwise R between lipids, 
nominally by penalising tail structure mismatch (which 
we treat as implicit in tail length mismatch 0 ) across 
the midplane. The particular mechanisms responsible for 
the direct coupling are not crucial to our model, however 
- for comparison with the literature we can simply esti¬ 
mate an effective strength of the conventional inter-leaflet 
mismatch energy 7 Q, which is shown on Fig. [3] as well 
The hydrophobic “indirect” coupling J 
promotes pairwise AR, by penalising mismatch in the bi¬ 
layer thickness profile. We also define J = 4J, which 
appears in the mean-field approximation of Eq. [T]used to 
derive the coarse-grained free energy, n can be related 
to the area compression modulus [ll, and penalises 
variation from species-dependent ideal length. Weaker n 
means the species can more easily adapt their tail length 
and structure to one another’s presence. The mismatch 
parameter Aq = iso — iuo is cast as a length, but repre¬ 
sents both tail length and structure mismatch; it couples 
to both the indirect and direct inter-leaflet couplings, J 
and B. Once fiducial values of the parameters are set, 
varying J alone approximates changing the mismatch in 
tail leng^ but not structure (e.g., adding carbons to one 
species Q ), while varying B alone approximates varying 
the mismatch in tail structure, e.g., unsaturation. We ar¬ 
bitrarily choose iso > iuo- The reference total thickness 
do = iso iuo is irrelevant in the absence of an external 
field acting on bilayer thickness. 
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B. Simulation protocol 

We simulate a = A/" bilayer where £ = 100 (script 
letters refer to the entire simulated bilayer, as opposed to 
a local bilayer patch). The Kinetic Monte Carlo scheme 
[H resembles Kawasaki (spin-exchange) dynamics, gov¬ 
erned by the Hamiltonian (Eq. [T]): lipids exchange be¬ 
tween neighbouring lattice sites within their leaflet, thus 
mimicking diffusive evolution after a quench from a high- 
temperature, randomised initial configuration. Hydrody¬ 
namics are not included, but we do not expect this to sig¬ 
nificantly alter the conclusions; we address this further in 
the Discussion. 

The overall leaflet compositions are conserved, since 
we do not consider flip-flop or exchange with the solvent, 
and are given by 


$t(b) = 

V 


( 2 ) 


where is the total number of S lipids in the top 

(bottom) leaflet. R and AR lattice sites can interconvert 
{SU -^US ^ SS + UU) to alter the “degree of registra¬ 
tion” (microscopic transbilayer symmetry) 


A = 


^fss + Muu 

M 


( 3 ) 


where Nss is the total number of SS lattice sites, etc. 
A can vary in the range 


+ - 11 < A < 1 - 


( 4 ) 


Although the overall leaflet compositions are conserved, 
A is not [ij. Hence, the bilayer can phase-separate into 
either locally symmetric or asymmetric modes while the 
conserved overall leaflet compositions may, as here, be 
fully symmetric between leaflets. 

Here we focus on symmetric, equimolar overall leaflet 
compositions ^^ = ^^ = 0.5. By Eq. lUthe degree of regis¬ 
tration A can vary from 0 (full pairwise antiregistration) 
to 1 (full pairwise registration). 


C. Parameters 

The lattice spacing is a ^ 0.8 nm. We use k = 
3a“^/cBT, corresponding to ka ~ 60/cBTnm“^, in the 
range for lipid bilayers at 300 K [2Q|-[^ . A fiducial value 
of the indirect coupling parameter is J rsj 2a-^kBT Q, 
while B leads to an effective value of 7 (Eig. [3] sec¬ 
ondary axis), for which existing estimates vary widely 
(7 - 0.01 -IksTum-^ 

We choose a mismatch parameter Aq = 2a ^ 1 . 6 nm, 
somewhat larger than the typical length mismatch in 
phospholipid mixtures, for which a (registered) phase 
thickness mismatch < 2 nm [2^, [23 would imply Aq < 
1 nm. Large Aq couples to both the indirect (J) and di¬ 
rect {B) couplings, increasing the energetic driving forces 


for both antiregistration and registration and making the 
competing phases clearer to interpret in simulation. The 
phenomenology is qualitatively similar upon reducing Aq 
[ 3 , with the caveat that Aq affects the effective value of 7 
arising from a given B (Eq. [8j). In Section Owe mention 
the dependence of the nucleation theory on Aq. 

We use V = 0.6 /cbT in the mean-field theory, above 
the mean-field threshold Vq = 0.5 /cbT for phase sepa¬ 
ration in the absence of other couplings, and use V = 
0.9 IvbT in simulation, where the corresponding threshold 
is = 0 . 88 /cbT due to fluctuations [25[. Although 

we cannot expect precise quantitative agreement between 
the mean-field theory and simulation, in this regime it ap¬ 
pears the exact value of V is not crucial to the kinetics; 
for example, a value V = 0.9 /cbT in the mean-field the¬ 
ory would yield a similar predicted landscape of relative 
R/AR growth rate Auj to that calculated in Eig. inii. 


III. REGISTERED AND ANTIREGISTERED 
PHASES 

We now discuss bilayer domain symmetry and asym¬ 
metry in terms of competing bilayer phase coexistences, 
focusing on those relevant to the present work and briefly 
considering those relevant to bilayers of asymmetric over¬ 
all leaflet compositions (^^ 7 ^ ^^)- Upon coarse-graining 
the microscopic model Eq. [T] (Eig. [2]:), a local bilayer 
patch is characterised by locally-averaged top and bot¬ 
tom leaflet compositions We derive a 

mean-field free-energy density /(0^,0^) [l|pwhich yields 
a phase diagram in (0^,0^) space (Eig. [j) [l[. The coex¬ 
isting bilayer phases determine the local order parameter 
in each leaflet, given by the projection of a given tie-line 
endpoint onto the (jA or axis. 


A. Definition of registration and antiregistration 

Some experiments, particularly with asymmetric over¬ 
all leaflet compositions, may use separate fiuorophores 
to image domain morphology in each leaflet [26l-[^. In¬ 
terpreting these in ( 0 ^, 0 ^) space, and relating them to 
experiments that do not image the separate leaflets , 
requires care. We define a registered (R) bilayer phase 
as one in which both leaflets are dominated by the same 
species. Hence, in our model, an R phase is dominated 
by SS or UU lattice sites (Eig. [ 1 ]) such that most lipids 
face one of their own species in the apposing leaflet. An 
antiregistered (AR) bilayer phase is one where the leaflets 
are dominated by opposite species. Hence, under our def¬ 
inition, “registration” is a property of a given homoge¬ 
neous patch of the bilayer, describing approximate local 
compositional symmetry between the leaflets. 

An alternative definition of registration is sometimes 
used, which we call “colocalised enrichment” [26|, [^ . 
This describes a bilayer in which the regions of largest 
top-leaflet composition (relative to the average in that 
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leaflet) spatially superimpose on the regions of largest 
bottom-leaflet composition (jf". Colocalised enrichment 
is therefore a property of domain morphology over the 
entire bilayer, not of an individual bilayer phase. It re¬ 
quires a tie-line of finite positive slope in (0^,0^) space, 
so that: i) both leaflets contain domains of larger and 
smaller than average and ii) the domains of large 

(jf belong to the same bilayer phase as the domains of 
large (ffi and are thus spatially colocalised with them. 

R-AR coexistence (Fig. [1] dotted line) can be accessed 
by a bilayer of asymmetric overall leaflet compositions. 
R-AR tie-lines have positive slope, and thus exhibit colo- 
calised enrichment, although the AR phase is highly 
asymmetric in composition. Hence, the colocalised en¬ 
riched domains in each leaflet reported in [13, [HI (where 
the leaflets were separately imaged) are consistent with 
either an R-R or R-AR tie-line. Quantitative compo¬ 
sition information would be required to unambiguously 
determine which was observed. 

In the present work, we consider a bilayer of symmet¬ 
ric overall leaflet compositions for which R-R and AR- 
AR tie-lines compete. In this case, “registration” versus 
“colocalised enrichment” are practically equivalent. R-R 
has positive tie-line slope (colocalised enrichment), and 
both bilayer phases are compositionally symmetric (reg¬ 
istered phases, under our definition). Conversely, AR-AR 
tie-lines are negatively-sloped so that enrichment in one 
leaflet colocalises with depletion in the other, and both 
bilayer phases are compositionally asymmetric (antireg¬ 
istered phases, under our definition). 

The literature is ambiguous. For example, full colo¬ 
calised enrichment (which can reflect either an R-R or 
R-AR tie-line) is described in as “registration”. In 
contrast, in [^, “registration” is also defined as the pres¬ 
ence of Lo (or Ld) on both sides of the bilayer, which is 
more like our definition of R phases versus AR phases. 
This ambiguity is not surprising because often, as here, 
the two definitions outlined above are similar. The two 
concepts can clash for asymmetric overall leaflet compo¬ 
sitions [HI where R-AR tie-lines can play a role. 

In asymmetric supported bilayers in the leaflets 
were imaged with separate fluorophores, revealing colo- 
calised domains in both leaflets. The same was found 
in highly asymmetric vesicles in [HI- On the face of it, 
these findings differ from some supported bilayers in [U 
which exhibited compositionally asymmetric regions in¬ 
terpreted as domains in only one leaflet. These were in¬ 
ferred from variations in total bilayer thickness without 
imaging the separate leaflets. Height mismatch smaller 
than a known value for R-R was measured, indicating 
either R-AR or AR-AR-R coexistence. R-AR was then 
inferred by detecting that saturated lipids were predom¬ 
inantly in the top leaflet. The system was interpreted as 
having domains in only the top leaflet H. 

However, it is probable that the asymmetric bilayers in 
both [H, Hsl and [2^ represent an R-AR tie-line. Thus, 
imaging the leaflets in [^ separatel y c ould have shown 
colocalised enrichment just as in [H,®- The top-leaflet 


gel domains in [U may have apposed regions in the bot¬ 
tom leaflet weakly more “gel-like” than the average in the 
bottom leaflet. Similarly, the colocalised enrichment in 
[H, IHl does not require that bottom [2^ (or inner [HI) 
leaflet domains are truly Lq - for R-AR they will, like the 
rest of the bottom/inner leaflet, be dominated by unsat¬ 
urated lipids. It only requires that they are weakly more 
“Lo-like” than the average in their leaflet, i.e., that the 
R-AR tie-line is tilted. Hence, the question of whether 
domains in one leaflet “induce domains in the other” be¬ 
comes the question of whether the given tie-line is tilted 
enough for both leaflets to exhibit deteetable domain for¬ 
mation when a different fluorophore is used in each. The 
degree of R-AR tie-line tilt will depend on the direct 
inter-leaflet coupling B and hence on molecular features. 


B. Kinetics and competing phase coexistences 

For most typical parameter choices, the R phases are 
lower in /(0^,0^) than are the AR, hence R-R tie-lines 
are equilibrium and AR-AR metastable [l|. As well as 
the coarse-grained bulk free energy /((/)^, 0^), the under¬ 
lying microscopic model yields free-energy costs for com¬ 
position and thickness gradients^ which describe the role 
of spatial structure and domain formation in the kinetic 
competition of metastable and equilibrium states. Linear 
stability analysis of the initial homogeneous state (Fig. [3]) 
predicts growth rates of competing R versus AR insta¬ 
bility modes and thus whether symmetric or asymmet¬ 
ric domains form first [l|. The interplay between gradi¬ 
ent and bulk free energies then governs the nucleation of 
symmetric domains (introduced in Section [V|), which de¬ 
termines the eventual fate of a bilayer that has initially 
become metastably AR-AR. The ability to derive both 
gradient and bulk free energies from defined microscopic 
interactions is a key advantage of the present model over 
a purely phenomenological approach. 

In Fig. O R-R separation is equilibrium (except for a 
very small region Bofi/k^T < 0.005). Below the “AR 
instability” line, AR-AR separation is also possible, but 
metastable. Below the “R instability” line, the initial 
state is not subject to R instability although R-R sepa¬ 
ration remains the equilibrium state. The red and blue 
colours show the difference in growth rates, Acj, between 
R and AR instability modes. For example, red signifies 
a faster-growing AR instability mode, so that AR-AR 
separation dominates initial demixing. 

For physical parameter ranges [l| , neither the R or AR 
mode is trivially dominant, so that moderate changes to 
lipid tail length mismatch (affecting the effective value of 
J) or tail structure mismatch (affecting B) can determine 
whether the initial instability after a quench leads to lo¬ 
cally symmetric (R) or locally asymmetric (AR) domains. 
In addition, the long-lived AR-AR states simulated in 
ilS provide prima faeie evidence that metastable trap¬ 
ping due to failure to nucleate R domains is possible for 
physical phospholipids. However, due to small Simula- 
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Q #1 < I » I# ir Red/blue colours and dashed/dotted lines: theoretical predictions. 

0 12 3 4 


Indirect (hydrophobic mismatch) 
coupling Ja^ / k^T 


Dots: degree of registration A at the end of the simulation. 
Circled numbers and bold boundaries mark kinetic classes. 


FIG. 3. Red/blue colours and dashed/dotted lines: mean-field theory parameter map showing relative growth rates of AR 
versus R instability modes for a bilayer comprising equimolar mixed leaflets, from linear stability analysis of initial demixing 
[I[, with Ao — 2a^ — 3a~^kBT, V = 0.6 /cbT. The equilibrium state is R-R in all cases. Below the “AR instability” line 

AR-AR coexistence is possible but is metastable. Below the “R instability” line the homogeneous state is not unstable to the 
R mode, although R-R separation is still the equilibrium state. Overlaid dots from simulation (with V = 0.9/cbT) show the 
average degree of registration A at the end of the simulation time, A = 0 (black, AR-AR) to A = 1 (white, R-R). We identify 
three kinetic classes discussed in the text (circled 1, 2, 3, bold line marks approximate boundaries). Illustrative simulation 
snapshot sequences for each class are shown (£ = 200 in the snapshots). Simulations are visualised with OVITO [5^. 


tion sizes in such studies, this could constitute a stable, 
not metastably trapped, state, an issue which we address 
in the Discussion. 

Finally we note that, contrary to [s^, AR-AR coex¬ 
istence does not require an exactly equimolar (or equal 
area fractions) mixture in each leaflet. Firstly, an overall 
composition away from = 0.5 can lie on one of a 

set of AR-AR tie-lines running parallel to the central one 
depicted in Fig. mi. Secondly, even if the overall com¬ 
position is outside any AR-AR tie-line, AR-AR-R coexis¬ 
tence can occur, in which the presence of some R phase is 
forced by the composition being far from muimolar, but 
the R phase coexists with two AR phases Q . In that sit¬ 
uation, every region of the bilayer except the R phase is, 
as for AR-AR, locally asymmetric. Like AR-AR, for typ¬ 
ical parameters, AR-AR-R is a metastable state, which 
could become trapped for strong hydrophobic mismatch 
or stabilised if domain size is limited (e.g., by simulation 
size). 


IV. KINETIC PHASE DIAGRAM 

Fig. [3] shows the results of simulations of phase- 
transition kinetics for varying indirect coupling J and 
direct coupling B, to model varying lipid tail length mis¬ 
match and structural mismatch. The overlaid greyscale 


dots signify whether registered or antiregistered domains 
dominate at the end of the simulated time (t = 10^ Monte 
Carlo Steps) by measuring the degree of registration A, 
with A ~ 0 (black) corresponding to AR-AR coexistence 
and A ^ 1 (white) to R-R. Each dot on Fig. [3] is an 
average of four independent trajectories. 

This simulated “kinetic phase diagram” agrees semi- 
quantitively with the theoretical linear stability analysis 
of initial demixing (red/blue on Fig. [3]): inside the blue re¬ 
gion, where R-R coexistence should be accessed directly, 
the simulation exhibits full registration. If instead the 
AR mode is fastest (red), we find two possibilities - the 
bilayer may reach R-R coexistence by nucleation out of 
the AR-AR state, or remain metastably trapped in AR- 
AR. Fig. 2^ shows the successful formation of registered 
nuclei, which grow from the boundaries of antiregistered 
domains. In Fig. IUd, a different random quench in the 
simulation leads to failure to nucleate despite unchanged 
parameters, illustrating the stochastic nature of the nu¬ 
cleation process. 

Having identified the expected three classes of kinet¬ 
ics, we next model the energetics of nucleating registered 
domains, to clarify the fates of the metastable AR-AR 
state (the second and third kinetic classes identified in 
Fig. [3|). 
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FIG. 4. (a) Snapshots of a trajectory in which two regis¬ 

tered nuclei grow to reach equilibrium, (b) Trajectory with 
the same parameters but a different random initial configura¬ 
tion, in which no nuclei survive. Parameters: Ao = 2 a, n — 
3a V = 0.9/cbT, J = 4a ^/cbT, B = 0.42a ^/cbT, 

C = 100 . 


V. NUCLEATION THEORY 


The metastability of AR-AR implies that arbitrarily 
small composition fluctuations decay. To reach R-R, reg¬ 
istered domains must be formed by composition fluc¬ 
tuations sufficient to generate a nucleus of a registered 
phase large enough for the area-dependent payoff in bulk 
free energy to outweigh the penalty for hydrophobic mis¬ 
match at the edge. We now study the energetics of this 
nucleation. For simplicity we consider the R nucleus to 
be compositionally uniform (thus ignoring mixing effects 
at the domain boundary), and use the response of the 
thickness profile to mismatch at the nucleus’ boundary 
to calculate the nucleation energetics. 

Fig. 0^ exhibits typical nucleation of registered do¬ 
mains. Our goal is to model the role of hydropho¬ 
bic mismatch around the edge of R nuclei in determin¬ 
ing whether they successfully grow (Fig. |4^) or decay 
(Fig. St)). For the reference (zero free energy) state, we 
take a registered bilayer at equilibrium with no transmid¬ 
plane mismatch and domains coarse enough for edges to 
make a vanishing energy contribution. We then assume 
that a dominant AR instability mode has led to an ini¬ 
tially antiregistered state, which incurs everywhere an 
energy cost, relative to the registered reference state, due 
to the direct coupling B. We then introduce a circular 
R nucleus of radius R, thus removing the direct coupling 
energy within the domain’s area but introducing a thick¬ 
ness mismatch around its perimeter (Fig. [5]). 

In the continuum limit of the lattice model, the free 


energy is 


Geont = 11 ^°^ “ 

(5) 

where a is the lattice spacing and the J term captures 
the thickness gradient penalty from the corresponding 
neighbour interaction (Eq. [T]) in the limit of small lat¬ 
tice spacing. We neglect contributions from the Ising- 
like interaction V. It is thought that, as hydrophobic 
mismatch increases, it becomes the dominant contribu¬ 
tion to the line tension of (registered) domains [2^, [sTI 
(this does not imply that thickness mismatch uniquely 
determines phase-transition temperature Isl). However, 
in Appendix [B1 we briefly discuss an upper bound for the 
influence of V on the following nucleation energetics. 

The leaflets’ compositions are reflected in their ideal 
thicknesses ^Q^(r). The actual thickness profiles 
adopted by the leaflets minimise the free energy Gcont- 
The actual thickness difference and total thickness pro¬ 
files are denoted A(r) = £^(r) — £^(r) and d(r) = £^(r) + 
£^(r). The calculations underlying the following results 
are detailed in Appendix 



A. Antiregistered background 

The nucleus appears within a phase-separated AR-AR 
background. In practice, R domains typically form at 
AR-AR domain boundaries (Fig. |4]). The effects of this 
on nucleation are briefly discussed in Appendix [HI but 
for the following calculations we are ignoring the Ising 
contribution V (which provides the only source of line 
tension at AR-AR boundaries). Hence, the pattern and 
orientation of the initial AR domains is unimportant, and 
we assume a spatially uniform state with S lipids in the 
top leaflet, U in the bottom. For simplicity we assume 
strong compositional segregation so that the phases are 
approximately pure, with ideal leaflet thicknesses given 
by 


^t, AR 
to 

nh, AR 
to 



Minimisation of Eq. [5] yields 

A AR _ ^^0 

K + 2B ' 


(6a) 

(6b) 


(7) 


for the difference in actual leaflet thicknesses. The species 
in each leaflet adapt their tail lengths (hence degree of tail 
ordering 0) to one another’s presence, balancing tail 
stretching with the direct coupling energy. Relaxing the 
assumption of strong compositional segregation would 
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FIG. 5. Schematic structure of a registered SS nucleus in¬ 
troduced to an antiregistered background. The leaflet thick¬ 
nesses inside and outside the nucleus are calculated from 
Eq. [151 The midplane (dashed) shifts discontinuously at 
r = R to maintain smooth outer contours (cf. [3, [13]). S 
(light) and U (dark) lipids illustrating composition are super¬ 
imposed. Primes indicate that leaflet thicknesses outside the 
nucleus differ from those inside. 


not qualitatively affect the physics. The physical content 
of Eq. [7| is simply that transbilayer-mismatched lipids in 
an AR phase cause a free-energy cost; in strong compo¬ 
sitional segregation this energy density is (see Eq. IA5p 

’ 2a^(K + 2B) ^ ' 

Since no mismatches in total ideal thickness exist, is 
uniform: 


= IgQ + IjjQ = do . (9) 


B. Introducing a registered nucleus 


We consider the energy change induced by a registered 
nucleus of the S species (Fig. |5j), with larger ideal thick¬ 
ness than its surroundings (the energetics are identical 
for a nucleus of the shorter species U): 


nt, nuc. 
^0 




(r) = iso , 

iso if r < R 
iuo if r > R 


( 10 a) 

( 10 b) 


The top and bottom leaflets now have the same com¬ 
position (ideal thickness) within the nucleus, removing 
the inter-leaflet mismatch within the region r < R and 
leading to an area-dependent free energy payoff via the 
B term of Eq. |5l The actual thickness difference profile 
is now 


•(r) = 


if r < i? 


kAo 


( 11 ) 


(a) R/a = 5 (b) = 1.41 



EIG. 6. (a) Galculated profile of total thickness d^^^'(r) rel¬ 
ative to the AR background, when the nucleus is introduced. 
The ideal thickness relative to the surroundings is Ao = 2 a. 
Arrow: decreasing decay length (b) d^^^'{r) relative to the 
AR background and normalised by its ideal value inside the 
nucleus, Aq. Arrow: increasing nucleus size. 


The discontinuity arises from the discontinuous compo¬ 
sition at r = R and the fact that, in contrast to the total 
thickness (i(r), no energy in the Hamiltonian penalises 
variation in the thickness difference A(r). 

If the composition interface were not sharp one would 
consider the coupling between composition and thickness 
gradients. Moreover, one could also allow for sliding into 
antiregistration at a registered domain boundary, which 
smears out thickness mismatch These would 

not qualitatively alter the fact that the registered nu¬ 
cleus’ boundary experiences an energy cost from thick¬ 
ness mismatch, but could help reduce it. 

Calculating the total thickness profile after the nucleus 
is introduced yields 


^nuc.(^) _ 


0 

<1 

[1 - f A-. (f) 

'*(i)] 

0 

<1 

[f/. (f) K„ 1 

(i)] 


if r < R , 

( 12 ) 

if r > R . 


Here, In and Kn are nth order modified Bessel functions 
of the first and second kind respectively. Their spatial 
dependence is controlled by a decay length 


e = 


2Ja?/h 


(13) 


which quantifies the competition between hydrophobic 
mismatch and stretching. Eq. [12] shows how the mis¬ 
match in ideal total thickness at r = R distorts the 
thickness profile both inside and outside the registered 
nucleus. The decay length ^ controls the lateral distance 
over which d^^^'{r) is perturbed (Fig. [ 6 ^). 

For large the R nucleus does not reach its ideal thick¬ 
ness even at its centre, r/a = 0. This is an important 
point; for small nuclei R ^ the effective thickness 
mismatch between the nucleus and its surroundings is 
less than its value Aq in the limit of large domain size 
(Fig. [ 6 ) 3 ). The thickness at the centre of the nucleus is 
given by 





(14) 


^ K, T 2H 


if r > R . 


d""^-(0)=do + Ao 
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Hence, small domains have smaller hydrophobic mis¬ 
match so experience a smaller effective line tension. This 
dependence of domain height on domain size could be 
observed by atomic force microscopy. Since J and n will 
typically be of the same order of magnitude [l| such that 
^ ^ a, it may only be significant for very small domains. 
However, increasing thickness mismatch of registered do¬ 
mains as they grow has been reported in molecular sim¬ 
ulation [33I, as predicted here. 

After introducing the nucleus, the individual thick¬ 
nesses nuc.^^^ ^b, nuc.^^^ bottom 

leafiets are given by 

^t, nuc.^^) ^ 1 A^^"-(r)) , (15a) 

^b, nuc.^^) ^ 1 ^ 

Assuming no empty space inside the bilayer, the leafiet 
thicknesses fully determine the outer contours of the bi¬ 
layer once the midplane position is specified. Noting that 
(thus ^^‘^•(r) and discontinuous 

at the nucleus boundary r = i?, we expect this discon¬ 
tinuity to occur in the midplane position to maintain 
smooth outer contours and minimise hydrophobic expo¬ 
sure . This can also be seen in molecular simulations 
containing AR domains [9[. The resultant bilayer struc¬ 
ture is shown in Fig.[5l Outside the R nucleus, the leafiets 
adapt their thicknesses (thus tail structure) to one an¬ 
other’s presence (Eq. [11]), i.e. the S species are shorter 
than iso and the U longer than iuo- Inside the nucleus, 
the registered region is able to achieve its ideal thickness 
only once the nucleus grows large enough. 


C. Energy of a registered nucleus 

Given the thickness profiles before and after the nu¬ 
cleus is introduced, the energy required to create a regis¬ 
tered nucleus of size R is calculated as (see Appendix 

■ <“) 

The negative term oc R‘^ arises from the removal of 
inter-leafiet tail structure mismatch over the domain’s 
area, and the positive term arises from hydrophobic mis¬ 
match. For large nucleus size R relative to the decay 
length the approximations In{x) ^ exp(x)/\/27rx and 
Kn{x) ^ exp(—x)A/7r/2x for large x give ^ 

(,/2R. The hydrophobic mismatch term then becomes 
overall linear in R, acting like a standard line tension 
r = AQb:^/8a^. Hence, for large nuclei or large stiff¬ 
ness b:, AG behaves as in 2D classical nucleation theory 
(CNT). For a«0.8nm, and Ao~a 

at T = 300 K, we estimate T ^ 2 pN between R and AR 
domains or T 8pN between R-R domains of differ¬ 
ent species, quite close to existing estimates for phospho¬ 
lipids [sj. Thus, although our simplified model does not 



FIG. 7. Energy AG for a registered nucleus of radius R, cal¬ 
culated from Eq. [161 with Ao = 2a and k = ?>a~^k^T. We 
vary the decay length ^ by varying J but hold the ratio of 
indirect and direct couplings constant at J/B — 2 {J/B — 8 
line in Fig. El). In 2D classical nucleation theory the ratio 
AG/R would yield a straight line of negative slope. Circles 
mark critical radii R . Inset: behaviour of AG without nor¬ 
malisation by nucleus size. 


capture all details of hydrophobic mismatch and the con¬ 
comitant bilayer deformation [s^, the associated energy 
scale is well captured. 

Some illustrative nucleation energy curves are shown 
in Fig. 0 For GNT in 2D, the ratio AG/R would be a 
straight line of negative slope. Deviation from a straight 
line indicates deviation from GNT. If ^ is small (weak 
hydrophobic mismatch compared to stiffness), the be¬ 
haviour is GNT-like over most of R; moreover the critical 
radius R (defined as the value where AG(R) is max¬ 
imised) occurs deep into the GNT-like regime. For larger 
^ (stronger hydrophobic length mismatch relative to stiff¬ 
ness) the critical radius no longer occurs in the GNT-like 
regime, so non-CNT effects arising from nonlocal thick¬ 
ness deformation infiuence nucleation. 


D. Nucleation barrier and critical radius 
landscapes 

Fig. □ (inset) implies that increasing J for a fixed ra¬ 
tio J/B and fixed n (equivalently, increasing ^) first in¬ 
creases, then reduces, the nucleation barrier. In the limit 
of large k one would expect increasing J and B to simply 
increase the energy scale, hence the nucleation barrier. 
But when J is comparable to the mechanism illus¬ 
trated in Fig. [6)3 takes over; increasing ^ = yJ~2Ja?^ 
reduces the effective thickness mismatch for small nuclei, 
thus reducing the hydrophobic penalty for nucleation and 
reducing the nucleation barrier. 
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Indirect (hydrophobic mismatch) 
coupling Jsf / k^T 

FIG. 8. Simulated kinetic phase diagram (notation as in 
Fig. E]), overlaid on the theoretical nucleation barrier for an 
R domain nucleating from an AR state (Ao = 2a, k = 
3 a~‘^kBT). The arrowed line shows J/B = 8, as used in 

Fig. El 



0 12 3 4 

Indirect (hydrophobic mismatch) 
coupling / k^T 

FIG. 9. As Fig. [8] but showing the calculated critical radius 
for R domain nucleation. 


The nucleation theory can be related to the simulations 
by plotting landscapes of the nucleation barrier AG{R*) 
(Fig. [8]) and critical radius R* (Fig. [9]). A small barrier 
clearly facilitates nucleation, although the role of the crit¬ 
ical radius is more subtle |35l-[^ . One important factor is 
that for smaller critical radii, regions randomly enriched 
in registration from the initial quench are more likely to 
be near or above the critical radius, so can grow more 
easily. 

The predicted energetics for nucleation are consistent 
with the classes of kinetics identified in Fig.[3l Where the 
R instability mode is fastest, nucleation energetics are ir¬ 
relevant since the equilibrium phases are formed immedi¬ 
ately. If the AR mode is fastest, a successful subsequent 
transition to equilibrium R-R coexistence generally oc¬ 
curs where the predicted nucleation barrier and critical 
radius are smaller. As we have used a relatively large 
value of the mismatch parameter Aq = 2 a, we briefly 


mention the dependence on it [38|. The critical radii 
are independent of Aq, because Aq represents both tail 
length and structure mismatch so affects both the indi¬ 
rect and direct coupling strengths. Thus the only change 
to Fig.[9]would be in the effective value of 7 shown on the 
secondary axis, which scales as Aq (Eq. [8j). On Fig. [HI 
in addition to the change in effective 7, the nucleation 
barrier AG(i?*) scales as Aq. 


VI. DISCUSSION 

In this paper we have studied phase-separation kinetics 
in a bilayer subject to competing symmetric and asym¬ 
metric instabilities. We have modelled the nucleation of 
symmetric (registered, R) phases where the initial spin- 
odal instability has led to antiregistered AR-AR coex¬ 
istence (local asymmetry everywhere). We made use of 
a microscopic lattice model that allows study of both 
phase equilibria and phase-transit ion kinetics, through 
coarse-grained theory and direct simulation. For realis¬ 
tic parameters, the competing symmetric and asymmet¬ 
ric instabilities can be comparable in strength, so that 
changes to molecule properties can tip the balance [l| . 

If the fastest-growing instability mode is R, spinodal 
decomposition directly into equilibrium R-R coexistence 
takes place (class 1 Fig. [3]). If the AR mode grows fastest 
then a metastable AR-AR state forms, and nucleation is 
required to reach equilibrium (class 2 Fig. |3|). Unsuc¬ 
cessful nucleation leaves the bilayer metastably trapped 
(class 3 Fig. [3]). 

The key parameters for the initial demixing and subse¬ 
quent nucleation energetics are the indirect inter-leafiet 
coupling J via hydrophobic mismatch, favouring antireg¬ 
istration, and the direct coupling B which favours reg¬ 
istration. Physically, the hydrophobic tail length mis¬ 
match or structure mismatch would affect the effective 
values of J or 5 respectively. Both J and B couple to 
the stiffness k, which determines how easily mixed lipids 
can change their tail length and structure to adapt to one 
another’s presence. Here we have focused on an equimo¬ 
lar mixture in both leaflets such that AR-AR coexistence 
competes with R-R; the kinetic considerations are quali¬ 
tatively similar for other overall compositions subject to 
competing instability modes, in which other states in¬ 
cluding three-phase coexistence (AR-AR-R or R-R-AR) 
can enter 0. 

Nucleation implies a critical radius, which a registered 
domain must exceed so that the penalty of thickness mis¬ 
match at the perimeter does not outweigh the bulk free 
energy gain. This does not mean domains automatically 
become roistered beyond a certain size, as is sometimes 
implied [Ef [13, [HI . If nucleation is energetically pro¬ 
hibitive, the mere presence of large coarsening antireg¬ 
istered domains does not guarantee registered domains 
will ever form, just as assembling a large volume of su¬ 
percooled water does not guarantee a supercritical ice 
nucleus will form. That said, it is intriguing to con- 
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sider whether some hydrodynamic or curvature-mediated 
mechanism (neither of which are included here) could 
dynamically assist very large antiregistered domains to 
overcome the nucleation barrier and become registered. 

In the present simulations it appears that a nucleation 
barrier AG{R*) > 5/cbT inhibits nucleation on the sim¬ 
ulated timescale for this simulation size. In principle the 
nucleation rate / (per unit area per unit time) is deter¬ 
mined by AG(i?*) via I = /q exp (—AG(i?*)//cBT), al¬ 
though the unknown kinetic prefactor Iq severely limits 
the quantitative utility of such a picture. We also expect 
the critical radius i?* to determine how large random re¬ 
gions of inter-leaflet symmetry during the initial quench 
must be in order for them to grow. The critical radius 
also captures a size-dependence of potential relevance to 
cell membrane rafts and clusters, discussed below. 

By examining the effect of thickness mismatch at the 
boundary of a nucleating registered domain, we showed 
that the process departs strongly from classical nucle¬ 
ation theory due to deformation of the thickness pro¬ 
file, which reduces the effective line tension of small nu¬ 
clei. For small nuclei, the deformation essentially spreads 
over the whole nucleus, reducing the effective hydropho¬ 
bic mismatch. To our knowledge this behaviour has not 
been noted in previous theories which tended to focus 
on domains larg e enough to treat the boundary as a 
straight line [^, . The prediction is supported by a 
recent molecular simulation of coarsening of a quenched 
bilayer [^, where registered liquid-ordered and liquid- 
disordered domains became respectively thicker and thin¬ 
ner as they grew through time. 

The simulation method does not include hydrodynam¬ 
ics, which are expected to dominate domain coarsen¬ 
ing beyond a lengthscale ^ 10“^ m such that the 
purely diffusive dynamics simulated here would no longer 
apply. However, this lengthscale is far beyond both 
the size simulated here and the predicted critical radii 
(Fig. El). Note that hydrodynamics cannot influence the 
free-energy landscape of the system, and thus cannot 
change the competing metastable and equilibrium states. 

The phase equilibria of “macroscopic” domains (which 
coarsen until limited by bilayer size) are not influenced 
by edge energies and, for most reasonable parameters, we 
predict that R phases are lower in free energy so that R-R 
is the equilibrium state in bu lk H , as seen in fluorescence 
studies of large domains [^, . At the opposite end of 

the size spectrum, experiments reveal pairwise antiregis¬ 
tration at the single-lipid level [Il|,|43, as also reported in 
simulation [8[ and predicted by our mean-field theory Q • 
In between these regimes, edge energies can influence the 
equilibrium state; for example, the AR-AR coexistence 
in 1,11^ is probably metastable in the limit of large size 
but could be stabilised if the simulation box is too small 
to accommodate a supercritical registered domain. 

This has important physical consequences, because cell 
membrane rafts or clusters are not macroscopic. Small 
domain size in vivo could be due to elastic repulsion Q , 


hybrid lipids [i5|, critical fluctuations [i6[, active recy¬ 
cling [ 13,111 or another mechanism. In either case, it is 
crucial to recognise that edge energies could influence the 
thermodynamic preference and even stabilise an AR-AR 
type of domain formation, in which one leaflet’s local en¬ 
richment in longer species (relative to the background in 
that leaflet) colocalises with a relative depletion of such 
species in the other leaflet. 

Although the cell membrane is maintained out of equi¬ 
librium, thermodynamic driving forces can be expected 
to play a role. The potential importance of this size- 
dependence is underlined by the fact that the estimated 
critical radii (Fig. [9]) can be of the order of putative 
lipid raft sizes [^. The basic biophysical question is: 
if a cluster of longer lipids and proteins exists in one 
leaflet, does it colocalise a similar cluster in the oppo¬ 
site leaflet to maintain transbilayer structural similarity, 
or does it choose shorter lipids and proteins to maintain 
uniform thickness? Our work implies that finite-size ef¬ 
fects, metastable states and phase-transition kinetics can 
be key in determining the answer. 


Appendix A: Calculations for nucleation theory 

In the nucleation theory the response of the actual pro¬ 
files ^^’^(r) to variations in the ideal thicknesses is 

calculated as follows. Varying Gcont (Eq. E]) with respect 
to and gives a pair of differential equations that 
can be combined to yield independent equations for the 
thickness difference A(r) = t{r) — £^{r) and the total 
d{r): 


n (A(r) — Ao(r)) + 2BA{r) = 0 , (Ala) 

VK, {dQ{r) — d{r)) -1- 2Ja^-^ ^ = 0 . (Alb) 

dr \ dr J 


In the above we have used the ideal thickness difference 
and total profiles, Ao(r) = '^o(^) ~ ^o(^) — 

(-oir) + 

In the AR background the S and U species’ lengths are 
not equal to their ideal values - the direct coupling B en¬ 
courages equality of tail length (thus degree of ordering) 
across the bilayer. Eq. lAlal is solved by 


A(r) = 


^Ag(r) ^ 
Hv T 2B 


(A2) 


Inserting Eq.[6]for the ideal thicknesses yields Eq.[71 The 
total thickness before the nucleus is introduced is uniform 

(Eq.©. 

After introducing the R nucleus, the difference in ideal 
leaflet thicknesses is zero within the R nucleus (AQ(r) = 0 
for r < i?), which gives Eq. [TTl The total thickness 
profile is found by solving Eq. lAlbl (which is express¬ 
ible as a modified Bessel’s equation of order zero). We 
require that the gradient vanishes at the centre of the 
nucleus {d d{r)/dr\r=o = 0), and that the total thick¬ 
ness approaches its ideal value away from the domain 
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{d{oo) = iso + iuo = do). The profile d{r) and its gra¬ 
dient are required to be continuous at r = i? in order to 
match the profiles inside and outside the nucleus. With 
ideal thicknesses given by Eq. [TOl this yields Eq. [121 
To calculate the energy AG required to introduce the 
registered nucleus, we consider the initial energy 
from inserting Eqs. 0 and [9] into Eq. O and the en¬ 
ergy when the nucleus is introduced, from inserting 

Eqs. Hi] and Ha instead. The energy for introducing the 
nucleus is then 


AG = G™=t - ■ (A3) 

Since Eqs. [TT] and [12] split at the nucleus boundary r = R, 
we evaluate separate contributions to AG outside and 
inside the nucleus and find 






+ J(d--(r)-dof , (A4) 




2 \dr 

+ 'iid’^'^^-{r)-2esof - 


K , .2 A^kB 


2{k + 2B) 


(A5) 


AG°^^ contains only positive contributions, resulting 
from deformation of the total thickness profile by the 
nucleus. AG^^ contains similar deformation terms, but 
also a negative term resulting from the fact that there is 
now no inter-leafiet mismatch in the region r < R. The 
sum AG^^ -j- AG°^^ leads to Eq. [T6l 

Note that the negative term in AG^^ yields the es¬ 
timated inter-leafiet mismatch energy density 7 under 
our microscopic definition (Eq. [8]) [l|. Near the strong- 
segregation regime this is a close approximation to the 
difference in bulk free-energy density between registered 
and antiregistered phases (which can be used to construct 
an alternative “macroscopic” definition of 7 [I|), becom¬ 
ing exact in the strong-segregation limit assumed in Sec¬ 
tion E However, in principle one could replace this term 
with the actual free-energy difference in cases where the 
correspondence between microscopic and macroscopic 7 
breaks down [I|. 


Appendix B: Effect of V on nucleation 

In Section m we ignored the effect of V (headgroup 
interactions) on the nucleation energetics. This applies 
best where hydrophobic mismatch is the dominant source 
of line tension j^. Here we sketch an upper bound on 
the contribution of V to nucleation energetics of a regis¬ 
tered domain. Eig.[TQl shows an R domain nucleating at a 


before nucleus 



with nucleus 

R 




FIG. 10. Schematic of an ASS' R nucleus of radius R forming 
at the boundary between AR-AR domains. 


boundary between AR domains (cf. Eig.[l|). This relieves 
some AR-AR interfacial energy cost (like a particle in a 
Pickering emulsion). 

Assuming the original AR-AR boundary is straight 
and the R nucleus is a circle, the length of AR-AR bound¬ 
ary removed by an R domain of radius R is equal to 2R. 
The length of R-AR boundary introduced by the nucleus 
is 2ttR. It is extremely difficult to estimate the line ten¬ 
sion caused by V, because i) the separated phases will 
not be fully pure and ii) the interface will relax its en¬ 
ergy by smearing the composition change over some finite 
distance, and iii) this will be coupled to the thickness 
gradients. However, an upper bound can be obtained 
by assuming, as in the main text, that the phases are 
strongly-segregated and that the compositional interface 
is sharp. With these assumptions, the line tension of the 
AR-AR boundary is 2E/a, and the contribution from V 
to an R-AR boundary is V/a. Subtracting the AR-AR 
energy contribution and adding the R-AR contribution, 
the correction to Eq. [16] would be: 

2V 

AG AG -f —(tt - 2)R . (Bl) 

Hence, the contribution of V is always to increase the 
nucleation energy (thus both the critical radius and nu¬ 
cleation barrier) of registered domains. Eor the param¬ 
eters ranges covered by Eig. [3] (assuming V = 0.9 /cbT 
as used in the simulations), we find that including the 
upper-bound correction from V increases the critical ra¬ 
dius R* by at most ^ 40% for J > 2a~‘^k^T^ from the 
corresponding value on Eig. [9] (If Aq is reduced to 1 a, 
reflecting decreased tail hydrophobic and structural mis¬ 
match so that V is proportionally more important rel¬ 
ative to the indirect and direct couplings, the increase 
can be ^ 100%). Eor AG(R*), if we assume that the 
original R* is not altered much by including E, the as¬ 
sociated increase in the nucleation barrier is estimated 
as (2V/a){'K — 2)R*. Alternatively, explicitly taking into 
account the change in R* resulting from including E, the 
increase in AG(R*) in the parameter range covered by 
Eig. [3] is up to ^ 100% for J > 2a~‘^k^T^ from the 
corresponding value in Eig. [8] 

We reiterate that Eq. IBII must overestimate the contri¬ 
bution E to nucleation energetics, possibly quite severely, 
so that the argument outlined here provides only a rough 
upper bound. 
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